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Abstract 

It is shown that there is an overlooked mechanism whereby some kinds of dissipation 
can enhance the Benjamin-Feir instability of water waves. This observation is new, and 
although it is counterintuitive, it is due to the fact that the Benjamin-Feir instability 
involves the collision of modes with opposite energy sign (relative to the carrier wave), 
and it is the negative energy perturbations which are enhanced. 



1 



The discovery of the Benjamin- Feir (BF) instability of travelling waves was a 
milestone in the history of water waves. Before 1960 the idea that a Stokes wave 
could be unstable does not appear to be given much thought. The possibility that 

n □ n fi 

the Stokes wave could be unstable was pointed out in the early 1960s [1|, 12, 131, |4j, but 
it was the seminal work of Benjamin and Feir which combined experimental 

evidence with a weakly nonlinear theory that convinced the scientific community. 

Indeed, Benjamin & Feir started their experiments in 1963 assuming that the 
Stokes wave was stable. After several frustrating years watching their waves disin- 
tegrate - in spite of equipment and laboratory changes and improvements - they 
finally came to the conclusion that they were witnessing a new kind of instability. 
The appearance of "sidebands" in the experiments suggested the form that the per- 
turbations should take. A history of these experiments and the outcome is reported 
in fl. 

The theory of the BF instability is based on inviscid fluid mechanics, and the 
assumption that the system is conservative. Therefore it is natural to study the 
implication of perturbations on the system. The implications of a range of pertur- 
bations on the BF instability have been studied in the literature: for example the 



effect of wind js, 9] and the effect of viscosity @, 10, 11, 12, kJ. Some perturbations 
have been shown to stabilize and others destabilize the BF instability. 

However, there is a fundamental overlooked mechanism in all this work. Mathe- 
matically, the BF instability can be characterized as a collision of two pairs of purely 
imaginary eigenvalues of opposite energy sign as shown in Figure m in y, this 
observation is implicit but the demonstration and implications have not been given 
heretofore. This characterization of the BF instability also appears in the nonlin- 



ear Schrodinger (NLS) model for modulation of dispersive travelling waves 15|, Il6 |. 
The eigenvalue with smaller positive imaginary part in the figure - just before col- 
lision - has negative energy, whereas the eigenvalue with larger imaginary part has 
positive energy. This energy is relative to the energy of the carrier wave E o]ies : 

E _ < ^Stokes < E+ where E± 

are the energies of the modes associated with the 
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FIG. 1: Schematic of the eigenvalue movement associated with the BF instability, as a 
function of the amplitude of the basic carrier wave (Amp) and the sideband wavenumber 
a. For a fixed a (vertical dashed line), there is a threshold amplitude. Below the threshold, 
the eigenvalues are purely imaginary. At the threshold a collision occurs and above the 
threshold, the eigenvalues are complex. 



respective purely imaginary eigenvalues in Figure [TJ Hence "negative energy" means 
that E_ - E Stokcs < 0. 

Once these facts are established we can appeal to the result that dissipation can 



destabilize negative energy modes 



13, 
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191 ] . There are many examples in fluid 



mechanics where negative energy modes - which are stable in the inviscid limit 
are destabilized by the addition of dissipation [^(J: Kelvin- Helmholtz instability 
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211 ] . interaction of a fluid with a flexible boundary 22j, |23J], stability of a fluid- 



loaded elastic plate [24| , Euler modes perturbed by the Navier-Stokes equations [25] . 

The book of Fabrikant & Stepanyants 21[ reports on experimental results for 
interfacial waves near the Kelvin-Helmholtz instability illustrating the enhancement 
of the instability of negative energy waves due to dissipation. See Figure 3.5 on page 
83, and the discussion on pp. 82-83 in [21]. 

In this brief communication, we sketch the basic result for water waves and then 
use a nonlinear Schrodinger model perturbed by dissipation for illustration. The 
NLS equation has shortcomings (e.g. symmetry which enables the phase to be fac- 



tored out, lack of validity for all time 



26j) but it provides a simple example of the 
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phenomenon 



Davey [10| gives a general argument for the form of a dissipation-perturbed NLS 
model, and Blennerhassett [8| starts with the full Navier-Stokes equations for a free- 
surface flow with viscous free-surface boundary conditions and derives a similar per- 
turbed NLS equation. For the dissipatively perturbed Stokes wave in deep water, 
these NLS models take the form 

\A t + (a - ia)A xx + \bA + (7 + ic) \A\ 2 A = , (1) 

where A is the envelope of the wave carrier, and the modulations are restricted to one 
space dimension x. When a = b = c = 0, equation (JTJ reduces to the NLS equation 
for the modulations of Stokes waves in deep water; hence a and 7 are positive real 
numbers. This NLS model has a BF instability and one can show explicitly that it 
involves a collision of eigenvalues of the form shown in Figure HJ We show below 
that when a > 0, there is always dissipation induced instability (before the BF 
instability), no matter how small a is. The parameter a is the perturbation of the 
rate of change of the group velocity dc g /dk due to dissipation. 

First, consider the linear stability problem for gravity waves in deep water. As 
the wave amplitude increases we show that there is a threshold value at which two 
eigenvalues of the linear stability problem collide, and these two modes have negative 
and positive energy. 

With 6 = x — ct, the speed c and amplitude r\ of the basic gravity wave of 
wavelength 2ir/k, to leading order, are c = co(l + k 2 e 2 + •••), Cq = gk~ l , 

V (9) =er h {9)+ e 2 m (6) + 0(e 3 ) , 

where e is a measure of the amplitude, 

m (9) = V2cos(k6 - O ) , 773(0) = k cos(2k6 - 26 ) , 

with 0o an arbitrary phase shift. Using standard results on integral properties of 
Stokes waves, the total energy relative to the moving frame is 

^stokes = T+v _ cI = v _ T ugi 2T = cl , 



where T and V are the kinetic and potential energies respectively, and / is the 
momentum 27]. Substitution of the Stokes expansion shows that _E stokes = o + 
e 3 Es + 0(e 4 ). Although the actual value of £3 is not important for the argument 
below, it is noteworthy that it is negative, and, using Table 2 of j^], one can confirm 
that E stokes is negative at finite amplitude as well. 

To formulate the linear stability for gravity waves take 

rj(9,x,t) i-> r}(9,e)+r](6,x,t) , 

where rf(6, e) is the basic carrier wave. Take a similar expression for the velocity 
potential <f>(9,x,y,t), where y denotes the vertical space dimension. Next one sub- 
stitutes this form into the water wave equations, linearizes about the carrier wave 
and takes rj(9, x, t) of the form 

T](9,x,t) = Re(E(#,cr)e kra ' +A *) , 

where a is real (the modulation wavenumber), and E(0, a) is periodic of the same 
period as the Stokes wave. The result is an eigenvalue problem for the eigenfunction 
E and eigenvalue A. 

The BF instability corresponds to a solution of this eigenvalue problem with 
< d < 1 and Re(A) > 0. When a is fixed - but nonzero and small - and the 
amplitude of the Stokes wave is increased, there is a threshhold amplitude where the 
BF instability occurs, and it corresponds to a collision of two eigenvalues as shown in 
Figured! To leading order the eigenvalues collide at A = ±ic 9 a, where c g = \\/g/k 
is the group velocity. 

To show that the colliding modes have opposite energy sign, we need a definition of 
the energy of the perturbation. This definition requires some consideration because 
the perturbation is quasiperiodic in space: 2ix/k— periodic in 9 and 2-njo— periodic 
in x. The total energy relative to the moving frame is 



i-2-n/a 1 p2w/k 

E tM =— / — / £dfldx. 
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where E = T — V — cl, V = ^gt] 2 , 



T= f + 20,0 X + + <%)6y , and I = f (0 e + (j) x )dy . 

J —oo J —oo 

Evaluating the perturbation energy for the two modes that collide leads to E total 

^Stokes + £ 2 E ± + . . . ? with 



E± = 2(k±a) (l-^l=F~) |C±| 2 . 

Here C± are scale factors associated with the eigenfunctions. Clearly sign(£' 2 f "£' 2 _ ) < 
for < a < 1. 

Having shown that the colliding modes have opposite energy sign, we consider a 
simple example which illustrates the mechanism for destabilization of negative energy 
modes by damping. A prototype for a conservative system, where the linearization 
has a collision of eigenvalues of opposite energy sign, which is perturbed by Rayleigh 
damping is 



q tt + 2bJq t + { X - r 2 )q + 25q t = , qG 



where r > is the "gyroscopic coefficient", \ a real parameter with \y\ r 2 , and 
5 > 0. 

The energy of the system (T5]) is strictly decreasing when 5 > and ||q t || > 0. Let 
q(i) = qe A '; then substitution into (j2J) leads to the roots 




A = ir - 5 ± Vx + 2ir5 - 5 2 and A = -ir - 5 ± Vx - 2ir5 - 5 2 . (3) 

When 5 = there are four roots A = ±i(r ± y/x)- The eigenvalue movement shown 
in Figure Q] is realized as x decreases from a positive value to a negative value, the 
collision occurring at % = 0. Suppose now x is small and positive (just before the 
collision) and look at the effect of dissipation on the two modes Ao = ir ± i^/x- 
Substitution of the eigenfunctions for these two eigenvalues into the energy shows 
that the mode associated with ir— iy/x has negative energy while the mode associated 
with ir + iy/x has positive energy. 

6 



FIG. 2: Schematic of the effect of dissipation on the eigenvalues associated with ([2]). 



With 5 small, expand the first pair of roots in ([3]) in a Taylor series 
X(5) = ir ± iy/x T -j=(t ± y/x) + <^ 2 ) • 

With < 5 <C 1 the eigenvalues are perturbed as shown to the right in Figure El 
The negative energy mode, Ao = i(r — y/x)> nas positive real part when dissipatively 
perturbed, and the positive energy mode, Ao = i(r + v/x), has negative real part 
under perturbation. Consequently, when small dissipation is added to the otherwise 
stable system (that is, < x ^ r2 )> the mode with negative energy will destabilize. 
After the collision (when % < 0) the growth rate of the instability is enhanced. 

It should be noted that other mathematically consistent forms of damping can be 
used. For example the uniform damping 

makes mathematical sense. But it leads to uniform contraction of the phase space, 
and does not destabilize negative energy modes. 

In order to study the effect of dissipation on water waves, one could start with 
the Navier-Stokes equations and perturb about the Stokes wave solution (see [s| for 
instance for the case of wind forcing). Another approach is to add viscous perturba- 



tions to the potential flow in various forms 



28] . From the modified equations one can 



derive a dissipative NLS equation. There are two issues to highlight: negative energy 
modes can be destabilized and so the BF instability can be enhanced by dissipation, 
and secondly, the form of the damping is important. It is known that negative energy 



modes of the Euler equations can be destabilized by the form of damping found in 
the Navier-Sto : 



Following 



[8 



£es equations [25]. 



10], a general perturbed NLS equation for various types of physical 
situations can be written in the form (CQ). The parameters a, b and c are taken to 
be non-negative. When they are positive, they represent dissipative perturbations, 
since the norm of the solution is strictly decreasing in time when a 2 + b 2 + c 2 > 0. 

When a = b = c = 0, the resulting NLS equation is a Hamiltonian partial 
differential equation; with A = u\ + iu 2 and u = (u\, u 2 ), 

Ju t = Vi/(u) + aJu xx — bJu — c||u|| 2 Ju , (5) 

where J was defined in (T5]) and 

H{u) = [ [!aK|| 2 -i 7 ||u|| 4 ] dx. (6) 
Jr 

Let 9(x,t) = kx — ujt + 6*o, and consider the basic travelling wave solution to flS]) 
when dissipation is neglected, 

(cos 9 — sin 9 \ 
( 7 ) 
sin 9 cos 9 I 

Then u , u, k satisfy — ui + ak 2 = 7||u || 2 . 

It is assumed that the Stokes travelling wave exists for a sufficiently long time 
before any dissipation can affect it: dissipation is taken to be a second order effect. 

Next we check the energetics of the BF stability problem in NLS. Linearize the 
partial differential equation (jSJ) with a = b = c = about the basic travelling wave 
([7j). Letting u(x, t) = 'Ro( x ,t){ a o J r Y ( x jt))j substituting into the conservative version 
of (jSJ), linearizing about u , and simplifying leads to 

Jv t + 2akJv x + av xx + 2^(u , v)u = , (8) 

where (•, •) is the standard scalar product on IR 2 . 

The class of solutions of interest are solutions which are periodic in x with 
wavenumber a. The parameter a represents the sideband. The BF instability will 
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be associated with the limit |er| <C 1. Therefore let 

oo 

v(x, t) = |v (t) + 2_,( v n{t) cosnax + w n (t) smnax) . 

71=1 

Neglecting the a— independent modes (superharmonic instability), the 
o— dependent modes decouple into 4— dimensional subspaces for each n, and 
satisfy 

Jv n + 2aknaJw n — a{na) 2 w n + 27u Uq v n = 

(9) 

Jw n — 2aknaJv n — a(na) 2 w n + 27u Uq w n = . 
When the amplitude ||u || = 0, it is easy to show that all eigenvalues of the above 
system (i.e. taking solutions of the form e At and computing A) are purely imaginary. 
Considering all other parameters fixed, and increasing ||u ||, we find that there is a 
critical amplitude where the n — 1 mode becomes unstable first through a collision 
of eigenvalues of opposite signature. 

To analyze this instability, take n = 1 and study the reduced four dimensional 
system 

Jvi + 2a/ccrJwi — aa 2 w\ + 27UnUn Vi = 

(10) 
Jwi — 2akaJ\i — aa 2 w\ + 27U0UQ Wi = . 

To determine the spectrum, let (vi,wi) = (q, p)e A *. Then (A, a) are determined by 
roots of 

A(A, cr) = A 4 + 2(p 2 + 4k 2 a 2 a 2 )X 2 + (p 2 - Ak 2 a 2 a 2 ) 2 , 

where p 2 = a 2 a 4 — 2a7||u || 2 cr 2 . Suppose p 2 > 0, then all four roots are purely 
imaginary (see Figure [1]) and given by 

A = \2aka ± ip and A = — \2aka ± ip . 

These modes are purely imaginary as long as p 2 > 0; equivalently || u 1| 2 < a 2 a 2 . 
Since «7 > 0, the instability threshold is achieved when the amplitude reaches 

\oicr\ 
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At this threshold, a collision of eigenvalues occurs at the points A = ±2ikaa; see 
Figure [I] for a schematic of this collision. 

It will be assumed henceforth that k 7^ 0. Then instability is through a collision 
of eigenvalues of opposite energy sign, which reproduces the instability mechanism 
for the full water-wave problem. 

Purely imaginary eigenvalues of a Hamiltonian system have a sig nature associated 



with them, and this signature is related to the sign of the energy [17|, |2j], |30j . Collision 
of eigenvalues of opposite signature is a necessary condition for the collision resulting 
in instability. 

It is straightforward to compute the signature of the modes in the NLS model. 
Suppose that the amplitude ||u || of the basic state is smaller than the critical value 
( TTTj) for instability. Then there are two pairs of purely imaginary eigenvalues, and 
they each have a signature. Let us concentrate on the eigenvalues on the positive 
imaginary axis 

A = iQ± with Q± = c g o ± p , c g = 2ak . (12) 

Then 

Sign(fi±) = i(q, Jq) + i(p, Jp) , 

where the inner product is real in order to make the conjugation explicit. One can also 
show that this signature has the same sign as the energy perturbation restricted to 
this mode. A straightforward calculation shows that Sign(f2±) = ±4 when ||uo|| = 0. 
Since p 2 decreases as the amplitude increases, the two modes will have opposite 
signature for all ||uo|| between ||uo|| = and the point of collision. 

Now consider the effect of the damping terms. Consider the reduced system (11 01) 
for the BF instability with the abc— damping terms included: 

Jvx + 2akaJwi — acr 2 vx + 27u Uq v x + T>\ = 
Jwi — 2afccxJv 1 — aa 2 w 1 + 27U Uq Wi + T> 2 = , 

with 

T>i = 2kaawi + aa 2 Jv x + bJvi + 2c(u , vx) Ju 

(14) 

T> 2 = —2kacrvi + aa 2 Jw 1 + bJwi + 2c(u , W]) Ju . 
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Now, let (vi(i), Wi(t)) = (vi, wi)e Ai . Then the eigenvalue problem for the stability 
exponent reduces to studying the roots of a determinant showing (with the help of 
Maple) that the two roots in the upper half plane are given by 

A± = 2ikaa - (b + aa 2 + c||u || 2 ) ± iVS , (15) 

with 

S = Aiaa 3 ka - c 2 ||u || 4 - Ak 2 a 2 a 2 - 4ia<r£;7||uo|| 2 - 2a 7( x 2 ||uo|| 2 + aV . 
When a = b = c = 0, these stability exponents reduce to 



A± = 2ikaa ± iy a 2 a A — 2a7cr 2 ||u || 2 . 

Now suppose these two eigenvalues are purely imaginary: the amplitude ||uo|| is 
below the critical value (TTTT) . To determine the leading order effect of dissipation, 
expand ( IToT) in a Taylor series with respect to a, b and c and take the real part 

-p / \ n (2^ h 2 2aka(aa 2 -7||u || 2 ) 

Re(A±) = -{aa +b + c\\u \\ ) t === +■■■ ■ (16) 

ya^cr 4 — 2a7<T^||uo|| 

For any a > there is an open region of parameter space where these two real parts 
have opposite sign since their product to leading order is 



Re(A_)Re(A+) = (aa 2 + b + c\\u \\ 2 y — ^— 9 ' ,/ + ■ ■ 



4a 2 k 2 a 2 (aa 2 — 7 ||u | |2 ' i:i 
a 2 a 4 — 2a7a 2 ||uo| 12 



For any a, b, c with a ^ there is an open set of values of ||u || where this expression 
is strictly negative, showing that Re(A_) and Re(A + ) perturb in opposite directions. 
In this parameter regime the dissipation perturbs the negative energy mode as shown 
schematically in Figure [2j 

It is clear that when only the 6— term is present all eigenvalues shift to the left. 
Therefore the b— term does not produce any enhancement of the instability, in agree- 

[ I 

ment with [12 ] . This damping is analogous to the uniform damping in (jlj). It is the 
a— term which leads to enhancement. However the NLS is a simplified model for 
water waves. 
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In summary, the fundamental observation is that BF instability is associated with 
a collision of eigenvalues of positive and negative energy, and there are physically 
realizable forms of damping which enhance this instability. It remains to be seen 
how this effect can be revealed in laboratory experiments, in numerical experiments 
based on the full water-wave equations, and in the open ocean. 
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Figure 1. Schematic of the eigenvalue movement associated with the BF instabil- 
ity, as a function of the amplitude of the basic carrier wave (Amp) and the sideband 
wavenumber a. For a fixed a (vertical dashed line), there is a threshold ampli- 
tude. Below the threshold, the eigenvalues are purely imaginary. At the threshold a 
collision occurs and above the threshold, the eigenvalues are complex. 

Figure 2. Schematic of the effect of dissipation on the eigenvalues associated with 
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